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Abstract 

Using a representation of multichannel quantum defect theory in terms of a quan- 
tum Poincare map for bound Rydberg molecules, we apply Jung's scattering map to 
derive a generalized quantum map, that includes the continuum. We show, that this 
representation not only simplifies the understanding of the method, but moreover 
produces considerable numerical advantages. Finally we show under what circum- 
stances the usual semi-classical approximations yield satisfactory results. In partic- 
ular we see that singularities that cause problems in semi-classics are irrelevant to 
the quantum map. 
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1 Introduction 



Electronic states of molecules are called Rydberg states, as opposed to valence 
states (whether covalent or ionically bound), when an outer electron moves far 
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away from the remaining ionic core. These states form electronic series which 
converge towards the ionization hmit of the molecule. 

The starting point of the quantum analysis of such states was the Quantum 
Defect Theory (see e.g. the review article by Seaton [1]), established first for 
atoms. It was shown that, due to the non zero spatial extension of the ionic 
core, the levels near the ionization limit follow the hydrogenic Rydberg law 
En = — Ry/ {n + dy, with only a constant (or nearly so) shift d of the principal 
quantum number n, entitled Quantum Defect. Quantum Defect Theory was 
extended to a Multichannel Quantum Defect Theory (MQDT), for the case 
that there are several series which converge to nearby states of the ion, and 
interact strongly. It was shown that MQDT gives a unified theory of bound 
states, autoionizing states and electron-ion scattering cross sections. This the- 
ory depends only on a small number of parameters, basically one quantum 
defect per interacting series. Practically, all is solved with matrices whose size 
is the number of series, while "brute force" methods would in principle try 
to diagonalize a matrix which contains an infinite number of levels for each 
series. 

This theory was extended to molecules by Fano [2,3]. There are always many 
interacting series corresponding to the rotational states of the ionic core. In- 
deed the slow velocity of the core rotation leads to a splitting of the rotational 
states of the core which is of the same order of magnitude as the sphtting 
between high lying electronic Rydberg states. The novelty were the implica- 
tions of the anisotropy of the core. The effect of this anisotropy on the ionic 
potential decays faster with distance r than the point charge 1/r Coulomb po- 
tential, at least as or Fano showed that the key point of the analysis 
is the existence of a cut off distance Tq. Below this distance the motion of the 
outer electron is tightly bound to the direction of the ionic core, above it the 
two become independent. Many detailed studies have followed on moderately 
excited Rydberg states of molecules, see e.g. reviews in refs. [4,5]. 

A novelty resulting from the studies of very high lying states of molecular Ry- 
dberg series was the experimental observation of "clear zones" in such inter- 
acting series. It was soon understood that this so called "stroboscopic effect" 
corresponds to resonances between the period of rotation of the core and the 
period of the outer electron orbit [6,7]. In order to study the relationship of 
this phenomenon with classical and quantum chaos, the classical limit of the 
MQDT theory was estabhshed [8]. A Poincare surface of section was intro- 
duced, whose coordinates correspond to parameters of the molecular system 
when the electron leaves the sphere of radius tq. The studies showed that the 
"stroboscopic effect" corresponds to a periodic reestablishment of a nearly 
integrable phase space structure at each resonance, while the classical phase 
space is completely chaotic in between. 
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The impetus for the present study was given by papers of Bogomolny where he 
introduces a semi-classical method for the quantization of a classical Poincare 
surface of section Map (PM) [9,10,11], and some other follow up papers, e.g. 
[12]. The original semi-classical method is not unique and is prone to singular- 
ity problems, thereby raising the question about its relationship with a true 
quantization. We followed exactly the opposite route, starting from a purely 
quantum method (MQDT), establishing its classical limit, and at last studying 
Poincare surfaces of sections. Our system is ideal for the study of this problem 
because we know in advance the correct result. Furthermore our system has an 
unusual and non trivial geometry, the phase space being a sphere. For bound 
systems, a preliminary report of the interpretation of MQDT as a Quantum 
Poincare Map (QPM) was provided in ref. [13]. However, the relationship with 
the theory of Bogomolny was only a formal analogy of the final quantum and 
semi-classical formulae used to compute the levels. In this paper we study and 
compare step by step classical, semi-classical and quantum descriptions of the 
same system. 

Another important property of the MQDT is, that it describes bound and 
ionized states in a unified way. Hence, we may use it to study the quantization 
of the so called Jung Scattering Map (JSM) for ionized states [14]. Indeed 
Jung proposed to send back ionizing trajectories onto the molecule in a well 
prescribed way, in order to obtain a of compactification of the phase space for 
above threshold states. The role of the JSM in this context, has been touched 
upon previously to elucidate certain properties of the ^'-matrix [15]. In the 
light of the reinterpretation of MQDT as a QPM for bound systems it seems 
very attractive to attempt a comprehensive description, which extends the 
advantages outlined in ref. [13] to include both scattering, bound states and 
resonances. 

The purpose of this paper is to present a classical map which combines the PM 
for the part of the phase space corresponding to negative energies of the excited 
electron with the JSM for the positive energy region, and then to interpret 
MQDT as the quantization of the resulting symplectic map. This will provide 
basically two advantages. On one hand a conceptual advantage, as the method 
and its approximation can be formulated in the context of a classical map 
and its subsequent quantization. In particular, we will see that phase space 
representations of wavefunctions for open channels quantize on trajectories of 
the JSM in the same way as bound states and ionization resonances quantize 
on trajectories of the PM. On the other hand on a more practical level, the 
dynamics is shown to be semi-separable [12,16] which permits a decomposition 
into separable steps, and suggests to choose a surface of section such that it 
separates these steps. Consequently, the QPM for half an evolution step of the 
system contains the relevant information, and its diagonalization will permit 
to find eigenvalues, eigenvectors, the level dynamics, etc in a simplified way. 
This will be particularly useful when dealing with near degenerate levels or 
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resonances which occur frequently and cause considerable difficulties in the 
conventional approach [13,17]. 

In the next section we shall develop the generahzed classical PM, which in- 
cludes the JSM for the particular conditions of a Rydberg molecule. In the 
third section we shall present the main result of this paper. MQDT for bound 
states and open systems will be formulated as a generalized QPM. In both 
these sections attention will be given to relaxing the common approxima- 
tion to keep the absolute value of the electron angular momentum fixed. In 
section four we shall discuss applications and show numerical examples that 
demonstrate the practical advantages of this method and may also give some 
new insight in situations of quantum chaotic scattering. These examples will 
exclusively be given within the approximation of fixed absolute value of the 
electron angular momentum, mainly to show nice two dimensional plots of 
surfaces of section. As the QPM was originally proposed in a semi-classical 
approximation [9,10,11], we shall in section 5 compare the approximate results 
of semi-classics to the ones of the "exact" QPM. We shall see that singularities 
causing problems in semi-classics are absent in the QPM. Therefore the QPM 
can be used to test regularization procedures for semi-classics. 



2 Generalized classical Poincare map for Rydberg molecules 

2. 1 Principle 

A Rydberg molecule is composed of an ionic positively charged molecular 
core and an outer electron. It is specifically given the name "Rydberg" when 
one considers states just below or above the ionization threshold, where the 
outer electron performs large excursions far away from the core. This situa- 
tion allows to define a distance tq which separates the potential felt by the 
outer electron into a long range Coulomb part, and a short range or collision 
part. In the long range part, the outer electron only feels the spherically sym- 
metric — 1/r Coulomb potential of the core (we use throughout the "atomic 
units" (a.u.) e — h — m — 1). lis motion is comparatively slow because one 
considers states near the ionization limit, but it can be either bound or un- 
bound. The electron moves along a Kepler orbit whose plane and axis are fixed 
in the laboratory frame. Meanwhile, the ionic core rotates freely around its 
angular momentum N. In the collision area the electron feels the lower sym- 
metry, faster decreasing, part of the potential of the core (dipolar (oc 1/r^) or 
quadrupolar (oc 1/r^) for a diatomic molecule) and it is strongly accelerated. 
This strong acceleration leads to a short range collision motion, which is ap- 
proximately the same for a total energy slightly below or above the ionization 
threshold, since the difference between the corresponding velocities is negligi- 
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Fig. 1. (a) The outer eleetron follows a Kepler orbit, fixed in space, for large 
distances r > ro and changes the Kepler orbit for r < ro by interaction with the 
short range part of the ionic core potential, (b) When the electron ionizes the Jung's 
recipe is to feed back the electron which escapes to infinity by doing a step backwards 
in time without collision with the core (i.e. a pure Kepler hyperbola with reverted 
velocity), and then re launch the electron towards the molecule along the resulting 
asymptotic trajectory for the next collision. 

ble compared to the velocity produced by the acceleration. This is basically 
the reason why in quantum mechanics MQDT treats bound and unbound 
states in a unified way. During the collision there is exchange both of angular 
momentum and energy between the electron and the core. The exchange of 
angular momentum implies a change of the plane of the Kepler orbit, since the 
angular momentum of the electron L is perpendicular to its plane (Fig. 1(a)). 
The exchange of energy leads to a variation of the eccentricity of the orbit, 
and may turn it into an ionizing trajectory. In this case the Jung's recipe [14] 
is to feed back the electron which escapes to infinity by doing a step back- 
wards in time without collision with the core (i.e. a pure Coulomb step with 
reverted velocity), and then re launch the electron towards the molecule along 
the resulting asymptotic trajectory (Fig. 1(b)). By doing this, one obtains 
trajectories which cross infinitely many times the r — ro sphere. Poincare sur- 
faces of section are described by a set of parameters of the molecular system 
at those instants where the electron crosses outwards the sphere. To define 
these parameters (several sets are possible) we first need to define precisely 
the reference frames and count the number of independent parameters for the 
determination of the dimension of this Poincare surface of section. 

2.2 Reference frames 

In this work, we will use three reference frames, the laboratory frame and 
two molecular reference frames, entitled "quantum" and "classical" reference 
frames. They are shown in Fig. 2. 
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Fig. 2. Reference frames for the molecular system, (a) Laboratory frame. Its Oz axis 
is selected along the total angular momentum J (b) "Quantum" molecular reference 
frame. Obtained from the laboratory frame with a rotation of Eulcr angles 95^/, 6m^ 
0. Its OZq axis is along the molecular axis M. (c) "Classical" molecular reference 
frame. It is obtained from the "quantum" by a rotation around M which brings the 
OXc axis along the angular momentum N of the core. This frame is physically the 
most significant, and all other figures of this paper are drawn with this convention. 
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(a) Laboratory frame. Due to the global rotational invariancc which entails 
conservation of the total angular momentum J, we can select, without 
loss of generality, its Oz axis along the total angular momentum J. The 
other two axis are arbitrary (but fixed) . In this frame the axis M of the 
core has (varying) polar angles 9m and 0m- 

(b) "Quantum" molecular reference frame. It is obtained by a rotation of the 
laboratory frame by the Euler angles 6m, ^m, (as common in quantum 
mechanics [18] we use the y convention for the Euler angles [19, Appendix 
B]). We entitled it "quantum" because it is the frame used in the quan- 
tum theory of molecules. The third Euler angle is arbitrary for a diatomic 
molecule, and is sometimes selected as 7r/2 [20] for normalization reasons, 
which do not concern us here. Therefore, we prefer to stick to the con- 
ventions of the first paper by Fano [2]. We restrict to the diatomic case 
for simplicity, a third Euler angle would be necessary for a polyatomic 
core. In this frame J hes in the "vertical" plane OXqYq. The angular 
momentum of the core N for a diatomic molecule is perpendicular to the 
core axis M, because the moment of inertia of the core around it is zero. 
It is thus in the "horizontal" plane OXqYq. 

(c) "Classical" molecular reference frame. This frame is obtained from the 
"quantum" frame by a rotation around OZq — OZq such that the new 
OXc axis is along N. As will become obvious below, from the physical 
point of view this frame is the most appropriate for the visualization of 
the results. As in previous publications [8,13,15,17,21] we will use it in 
all figures of this paper. 

Formulas to go from one reference frame to the others are obtained from simple 
geometric arguments based on the previous remarks, and are established in 
Appendix A. 



2.3 Count of parameters and the dimension of the Poincare surface of section 

The configuration space is described by five coordinates in the laboratory 
frame. Two angles 9m, describe the position of the axis of the diatomic 
ionic core (a third would be necessary for a polyatomic core). Two angles 
Of., ^Pe and a radial distance r describe the position of the electron. Thus con- 
figuration space has a dimension of five and phase space has a dimension of 
ten. Conservation of the total angular momentum J reduces the dimension of 
phase space by 4, two for the Poisson commuting actions and J^, and two 
for the associated angles. The first angle can be selected as the polar angle 
ipj associated with J^, that is the angle of the projection of J onto the lab- 
oratory plane Oxy, which is constant (with our choice of the laboratory axis 
it is in fact undetermined but this is of no importance). The second can be 
a global rotation angle around J associated to J^, which is cyclic (but not 
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constant): see e.g. the analogous reduction to the {ji, j2, j3,Wi,W2,w^} action 
angle system for the pure Kepler problem in [19, Chap. 10-7]. Consequently, 
the dimension of the reduced phase space in the general case is six, and a 
Poincare surface of section has dimension four. 

The introduction of the auxiliary approximation that the modulus L of the 
angular momentum L of the electron is conserved during the collision re- 
duces the dimension of the phase space by two (L^ and the associated angle), 
and thus the Poincare surface of section has dimension two, which enables 
graphical representations. This approximation was used in our previous works 
[8,13,15,17,21], and is valid in the experimental works which motivated the 
theoretical study [6,7]. The reason for its validity is that the splitting between 
different electronic levels produced by the spherically symmetric part of the 
potential of the core is much larger than the splitting between sub levels of 
L^, produced by its anisotropic part. In the experiments, L takes low values, 
restricted to 0. .. 2 by a two steps laser excitation. Here, however we want to 
study the semi classical limit, where the angular momenta take large values. 
There are two possibilties to increase angular momenta. First we can increase 
them without changing other molecular parameters. This is experimentally 
possible since the application of R.F. fields [22] or of combined electric and 
magnetic fields [23] in atomic Rydberg states enable to climb the L ladder 
up to circular orbits. In this case the constant L approximation breaks down 
for higher L. A study of the consequences for this case has been performed in 
ref. [24]. However, this is not what we aim at, since it implies a complete change 
of the physics of the problem at hand. We want to gain insight into the physics 
of the experimental works in refs. [6,7], by establishing its semi-classical limit. 
Mathematically this amounts to letting — * while keeping all classical pa- 
rameters constant. Practically, since h is a constant of nature (taken as 1 in 
atomic units) , this requires the increase of all classical parameters while keep- 
ing constant their dimensionless ratios. For example, decreasing hhy a factor 
of m implies to multiplication of all momenta by m, angular L, J, N, and 
principal, n. The Kepler period Tg of the orbital motion of the electron, which 
is proportional to n"^ is thus multiplied by m^, the period T/v of the rotational 
motion of the core which is proportional to I/N (where / is the moment of 
inertia of the core) must be multiplied by the same factor, so that / scales as 
m^. Similarly, one must scale the ratio between the geometric size of the core 
and the De Broglie wavelength, keeping constant the core's shape parameters 
(ratio anisotropic / isotropic). This maintains constant the ratio of the clas- 
sical periods associated with L, due to the isotropic part of the short range 
potential, and with A = L • M, due to the torque produced by the anisotropic 
part of this short range potential. This maintains the validity of the approx- 
imation of constant L. Of course this needs proportionally larger molecules, 
more difficult experimentally to deal with, but not completely out of reach 
since recent experimental studies have found Rydberg states on molecules as 
large as Benzene Argon complexes [25], Diazacyclooctane (DABCO) and bis 
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(benzene) chromium (BBC) [26]. 

Thus, for practical reasons, we shaU continue to use the constant L approxi- 
mation in many formulae and all pictures presented in this article. Yet we shall 
show, especially in the present and the next sections, that these restrictions 
are not essential for the interpretation of MQDT as a QPM, and how formulae 
can be generalized. 

Finally notice that in ref. [21], we started with only seven parameters. The 
three missing parameters were the three angles associated to the three con- 
served momenta J^, and L^. This is legitimate in classical mechanics since 
they can be computed afterwards by a mere quadrature. For the compari- 
son between classical and quantum mechanics, which is the main purpose of 
this work, this is natural since these angles are totally undetermined due to 
Heisenberg's uncertainty principle for well defined values of the corresponding 
momenta. 



2.4 Parameterizations of the Poincare surface of section 

We first discuss the case of constant L. For a proper parameterization of a 
two dimensional Poincare surface of section, we need two conjugated param- 
eters, a momentum and an angle. We will not use the most obvious angles 
9e, (pe of the electron and their conjugated momenta, but angles associated 
with angular momenta Lzq of the electron and N of the core because in the 
quantum MQDT description of the molecule, these momenta enter naturally 
as (approximate) quantum numbers for the limiting cases [2] . There arc differ- 
ent choices, depending on whether we work in the molecular or the laboratory 
frame. 

(1) "Quantum" molecular frame. We use the projection A = Lzq = L ■ M 
of the electron angular momentum onto the core axis, and the angle 
(fi^y the azimuthal angle of the projection of L into the horizontal plane 
OXqYq, conjugated like the ji,wi pair in [19, eq. 10-139]. This choice 
of parameters implies that the Poincare surface of section is a sphere 
on which every point is defined by 9^ — arccos(A/L) and (^2, i.e. the 
direction of L in the molecular frame. 

(2) Laboratory frame. We use the modulus N of the angular momentum of 
the core (not one of its projection onto laboratory frame). Intuitively, 
the associated angle is the angle of rotation ipj^i of the molecular core M 
around the direction of N, and the only possible reference for this angle is 
given by the direction of the total angular momentum J. A proof of this 
intuition, namely a precise definition of this angle, formulae relating the 
pairs { Lzq, } and { N, ip^ }, and the proof that this transformation 
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is canonical is given in Appendix A.l. 
(3) "Classical" molecular frame. We use Lz^ = Lzq (the two Z axis coincide), 
and (fl, the angle referred to N in the "horizontal" plane (Fig. 2(c)). No- 
tice that this pair is not conjugated, as discussed in section 5.1. The 
consequence is that the Poincare map does not conserve the area of the 
sphere with these coordinates, a small effect for the choice of parameters 
we have used in this paper. We stick nevertheless to this choice for graph- 
ical representations, because it gives more intuitive physical pictures as 
shown below. 

Notice that with this choice of parameters, the direction of the plane of the 
Kepler orbit is well determined in the molecular frame, since it is perpendicular 
to L, and its eccentricity is well determined by L and the electron energy Eg, 
but that the direction of its main axis in this plane, i.e. the Laplace-Runge- 
Lenz vector A [19], is not. As explained in ref. [8], the reason for this is that 
no component of A , which defines this major axis, commutes with which 
is kept fixed. Our counting of parameters implies that the direction of A is 
not an independent parameter. 

Introducing the approximation that L is a constant has as a consequence that 
the information about the main axis direction is hidden in the angle associated 
to L^. Then, the choice of action angle variables corresponds to that typically 
used in the pure bound Kepler problem in celestial mechanics [19]. The actions 
arc ji = Lz, j2 — L and js = ^J—1/2E^ (i.e. "principal quantum number"). 
The associated angles are the "longitude of the ascending node" Q {— ipL+'^f'^ 
with our notations), the "argument of the perihelion" u, which is the direction 
of the main axis in the orbital plane (or the direction of A in this plane) we are 
looking for, and the "mean anomaly" which is proportional to time. Thus, the 
angle of the main axis of the ellipse, which is constant for pure Kepler motion, 
i.e. for Rydberg molecules when the electron is far away from the core, but 
varies during the collision, is part of the map only, if the approximation 
constant is omitted. 



2. 5 Description of the classical Poincare map 

We now briefly describe the PM, both in the laboratory and in the molecular 
frame. The dynamics is composed of two independent consecutive steps: 

(i) collision step: when r <S Tq the electron feels the short range anisotropic 
part of the potential. We suppose that this step is very short, so that the 
core axis remains flxed (impulse approximation). It is thus best described 
in the molecular frame. Due to the cylindrical symmetry of this potential 
the projection A = L • M of L onto the core axis M is conserved. With 
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the auxiliary approximation that Li^ is constant, this amounts to the 
cohision being described by a precession of L around M by an angle 
S^-Pl- Furthermore the invariance of the potential under reflection with 
respect to any plane containing M implies that 5</?l is odd with respect 
to 9l — > n — 9l- The simplest choice is 6ipL — KcosOl, which defines 
the strength K of the interaction. Notice that this 6l dependant rotation 
about the OZ molecular axis is a twist, not a global rotation of the sphere, 
implying that the dynamics generated by the map can be chaotic [27,28]. 
Since L has changed, at least in direction if we make the approximation 
constant, due to the the conservation of total angular momentum 
J = L + N, N changes both in direction and in magnitude. This has far 
reaching consequences. The change of the magnitude entails a change of 
the rotational energy of the core Ej^ = BN^ {2B is the reciprocal of the 
core moment of inertia / according to usual spectroscopic notation [29]). 
Consequently, due to the conservation of the total energy E the electron's 
energy Eg = E — En also changes. If Ee becomes positive the molecule 
ionizes, if it remains negative the resulting changes of Tg and T/v will 
change the parameters of the next Coulomb step. While the "quantum" 
molecular reference frame is fixed during the collision, in the "classical" 
molecular reference frame, its OXq axis, which is chosen along N, will 
change (Fig. 2(c)). This produces an additional phaseshift Sip'^^, which 
was called frame recoil in ref. [8] where explicit formulae are given. It 
can also be computed from eq. (A. 3) of Appendix A.l, which gives the 
angular position of N in the fixed "quantum" molecular frame, by taking 
the difference of positions before and after the Scpi kick, 
(ii) Coulomb (or free rotation) step: when r S> ro the electron only feels 
the spherically symmetric — 1/r part of the potential. It moves on a Ke- 
pler orbit. Its angular momentum L and its Laplacc-Rungc-Lcnz vector 
A are constant in the laboratory Oxyz reference frame. Meanwhile the 
molecular core rotates at constant speed around its angular momentum 
N which is perpendicular to its axis M. When seen in the "classical" 
molecular frame, with OZ along M and OXc along N (Fig. 2(c)), L 
(and A) seem to rotate the opposite way (clockwise) around OXc- The 
total angle of rotation is 5/3 = —27iTg/TN, i.e. the ratio of the periods of 
the electron orbit and of the core rotation. We have Tg = 27r(— 2£'e)^/^ 
a.u. and Tn = 2n/2BN a.u.. 

The crucial point of this model is the exchange of energy between the electron 
and the core during the collision step. If we neglected it the model would coin- 
cide with with a "kicked spin" model [30,31], which cannot lead to ionization. 

For the free rotation step we prefer to visualize the results in the "classical" 
molecular frame. Indeed in the "quantum" molecular frame the free rotation 
is about a direction of N in the molecular OXqYq plane (see Fig. 2(b)), which 
varies from step to step. We will visualize ionization on the Poincare sphere in 
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the "classical" frame by noting that positive electron energies Eg correspond 
to points where the total energy E is greater than the rotational energy E^. 
Using = (L + N)^, this leads to the inequality 

En = BN^ = B (^-L cos a + ^Jj^- sin^ < E (1) 

This condition only depends on the angle a between L and N and accordingly 

the positive electron energy region corresponds to a cap around the positive 
N axis, which is fixed in this frame (Fig. 2(c)). If we relax the condition 
constant, we have a four dimensional PM and formulae for the collision and 
the recoil are modified, but eq. (1) remains valid. 

Practically, there are two possibilities to compute the classical PM. In all our 
previous works, we considered it entirely in the "classical" molecular reference 
frame, by computing in turn free rotation, collision Sipi,, frame recoil, new 
values of N, En, Tn, E^, Te and plotting the results in the "classical" reference 
frame by using equations given in Appendix A.l. In the present work we 
alternate between laboratory frame and "quantum" reference frame, in a way 
summarized in the first three columns of table 1 (the last two will be filled 
below in sections 3 and 5). We choose this more complex, and thus slower, 
procedure, because as outlined in section 5, the coordinates in the "classical" 
reference frame are not canonically conjugate. This is of no harm for purely 
classical computations, but forbids the semi-classical analysis. Of course we 
have checked that both representations give the same results. 

Finally we will describe how to construct the combined Poincare-Jung scatter- 
ing map. The JSM is defined such that the electron which escapes to infinity 
is fed back by doing a step backwards in time without collision with the core 
(i.e. a pure Coulomb step with reverted velocity), and then re launching the 
electron towards the molecule along the resulting asymptotic trajectory for 
the next collision (Fig. 1(b)). Notice that L and A are constant in the lab- 
oratory frame during such a pure Coulomb step. For the constant case 
this amounts to starting from the same point on the Poincare sphere in the 
molecular frame ii Eg > after a collision, i.e. to skip the three first steps 
in Table 1. Indeed the rotation of the core is neglected during the impulse 
collision step, and during journeys to and from infinity of the electron these 
rotations cancel since they are made for r > tq, where the short range poten- 
tial is negligible. In the general case the "argument of the perihelion" uj is also 
conserved. The proper JSM is a map from ionized to ionized states, the PM 
a map from bound to bound states, but there is obviously no technical prob- 
lem to iterate a combined map from unbound to bound states and vice versa 
by selecting the JSM recipe to feed back ionizing trajectories. This seemingly 
strange feed back was suggested and discussed in ref. [14]. We shall show that 
it is indeed a good representation of the results in quantum mechanics under 
semi-classical conditions. 
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Table 1 

Poincare Map: summary. Molecular coordinates are in the "quantum" reference 
frame. 



Coords. 


Action 


Classical 


MQDT 


Semi Classical 
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eqs.(A.8,A.ll) 
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eqs.(A.14,A.15) 
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2K cos Ol 


g2i7r^A 


FsiLz, <Pl) = Lz^PL + ^.TTHLz 


Lz^L 













3 MQDT as a generalized quantum Poincare map 



3.1 The usual MQDT theory 

MQDT for molecules is a well established, time honored theory [2,3,4,5]. We 
will outline here only what is needed for our purposes. The theory closely 
parallels the preceding classical ideas, which originated from the classical limit 
of MQDT, but which were imphcitly present in Fano's seminal paper [2], albeit 
with a fully quantum way of writing down the theory. Wave functions have 
different forms for r < tq and r > tq and they are matched at r = tq. Going 
from r = to r = cxD we have 

(i) Collision step: r < tq. The electron's motion is rigidly coupled to the 
molecular axis M. Its speed exceeds by far the speed of the motions 
of the nuclei of the molecular ion and the Born-Oppcnhcimcr factoriza- 
tion of the molecular wave function applies (Hund's case (b) coupling 
type according to molecular spectroscopy nomenclature [29]). Due to the 
cylindrical symmetry of the potential the projection A = L • M of the 
electron's angular momentum onto the molecular axis M is a constant of 
the motion, whereas the squared angular momentum of the molecular 
ion rotation does not have a definite value. We may or may not have a 
definite value of depending on the use of the auxiliary approximation 
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= const. The angular part of the wave function is thus an eigen- 
function of A (in addition to J and Mj) denoted by x\^j^\9[,./^'^^Wi) 
[2], where (^g, ip'^ are the angular coordinates of the Rydberg electron in 
the "quantum" molecular reference frame OXqYqZq defined in Fig. 2(b) 
and Appendix A. For r ~ ro, i.e. outside of the molecular core, where 
the potential is approximately pure Coulomb, the radial part of the Ry- 
dberg electron wave function is a linear combination of two independent 
Coulomb functions. We choose the one regular at r = 0, s{Ee,r), as well 
as one irregular at r = 0, c{Ee,r), suitably selected for its asymptotic 
properties when r ^ oo [1] . The normalized total wave function reads [2] 
as 

^49'„ (^;, M; r) = X^/'^V^, V^'e, M) {s{E,, r) cos(7r//A) + c{Ee, r) sin(7r/XA)) . 

(2) 

This implicitely defines the "quantum defects" /^a- Notice that Ee de- 
pends on A. If the subsidiary approximation constant is not used, the 
X's are not eigenstates of L^, but can be developed in such a basis, im- 
plying a summation over L and an extra label L on the /x's and the E^s. 
The quantum defects can in principle be computed by ab initio electronic 
molecular calculations between r = and r = vq, but we will take them 
here as semi empirical parameters, as is frequently done. We will assume 
that they are independent of energy, in a restricted range above or below 
zero electron energy. This is justified by the fact that in the inner region 
the change of the total energy is negligible as compared to the electron 
kinetic energy in the Coulomb well. The wave function in the collision 
region should not depend significantly on such small variations of the 
total energy. For large distances there is a difference, leading to bound 
or unbound states. We will neglect refinements used in actual molecular 
calculations to compute levels with somewhat lower principal quantum 
numbers, for which one allows quantum defects to vary slowly with E^. 
(ii) Coulomb step: r > Tq. The electron is coupled to the isotropic part of the 
Coulomb potential only, and cannot exchange angular momentum with 
the core. The core angular momentum N and the electron angular mo- 
mentum L are thus separately conserved, as is J = L + N. The angular 
part of the wave function ^^^'^'^\6e, ^Pe,^) corresponds to the coupling 
by a Clebsch-Gordan coefficient of an electron wave function labeled by L 
and of a molecular core wave function labeled by N (Hund's case (d) cou- 
pling type according to molecular spectroscopy nomenclature[29]). Here 
{6e, ^Pe) arc the electron coordinates in the laboratory reference frame 
Oxyz. The radial part of the wave function of the Rydberg electron is 
again a linear combination of regular and irregular Coulomb wave func- 
tions. The total wave function reads 

(^e, fe, M; r) = ^g/'^^Ve, ¥>e, M) {s{E„ r) + c{E„ r) (In) , (3) 
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where Ee = E — BN{N+1) can be positive or negative and cjv and djv are 
to be determined from asymptotic conditions. Here again an additional 
summation over L is necessary when omitting the auxihary approxima- 
tion that be conserved. 

For r ~ ro we have one and the same wave function developed into two 
different bases (2) and (3). The electron angular parts correspond to a change 
from the molecular frame to the laboratory frame. The orthogonal matrix U 
which performs this change [2], 

N 

is proportional to a Clebsch-Gordan coefficient: 

Una -{L- AJA|LJAro)(-l)^-^+^(2 - Sj,o)K (5) 

if we use the definition of rotational wave functions given in Appendix A. 

We transform a general wave function = J2a ^a^a from the molecular to 
the laboratory frame by combining (2) and (4). The matching of the wave 
functions (2) and (3) yields 

*(^e, ^e, M; r) = ^ ^S/'^^Ve, <^e, M)vl/^(r) 

N 

*Ar(r) = Yl Una (s(Ee, r) cos(7r//A) + c{Ee, r) sin(7r//A)) (6) 

A 

taking into account that the implicit E^. dependance of the radial wave func- 
tions on A in eq. (2) and on in eq. (3) is negligible near r ~ Tq. This would 
be true neither for r -C Tq nor for r 3> Tq. 

The coefficients Ai^ in eq. (6) are obtained from the asymptotic properties of 
the wave function for r — > oo. These properties in turn depend on whether the 
electronic energy in each asymptotic channel labeled by A^, namely E^ = E — 
BN{N + 1), is positive (open channels) or negative (closed channels). In any 
case we define a principal quantum number (non integer) by E^ — — l/2i/|f. 
For closed channels i/jv is real, for open channels it is imaginary. 

In a closed channel, for the motion to be bounded, the exponentially increasing 
part of the Coulomb wave function must be zero, which implies that the radial 
part ^n{t) of the wave function in eq. (6) must be proportional to a suitable 
combination of regular and irregular Coulomb wave functions which is regular 
when r — > oo, i.e. of the form: [1,32]: 

*Ar(r) oc - cos(7rz/Ar) s{Ee, r) + sin(7rz/Ar) c{Ee, r), (7) 
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Combining equations (6) and (7) yields 



J2UNAsm{7r{uN + fiA))AA = (8) 

A 



In an open channel the Coulomb wave function may be decomposed into in- 
coming and outgoing Coulomb waves </7^ which have the asymptotic behavior 

[1]: 

93^(r) = -j= (c ± is) ~ (7rA;)^2 exp (±iC) for r — > oo (9) 

v2 

IX i 
(^kr- -Ltt + - \n(2kr) + arg(r(L + 1 - ik)), k = —. 

Scattering is defined with respect to pure Coulomb scattering, that is, the 
eigenphases r are measured relative to the Coulomb phases Q. The wave func- 
tion will have the correct asymptotic behavior in the limit r — > oo for the 
asymptotic channel if its radial part is of the form: 

*jv(r) oc ^Cn' {(p~{i'N;r)SNN' - V^^{i'N;r)SNN') , (10) 

N' 

where the sum runs only over all open channels, and the Cn' are constants to 
be related to the Aj^. This defines the scattering matrix S. It is interesting to 
consider the matrix T which diagonalizes S, i.e. 



Snn' = X! exp(2i7rT^) T^n' (H) 



Using eq. (9), we obtain: 
I 

\ ^ N' J 

f^e {s{Ee, r) cos(7rT^) + c{E^, r) sin(7rT^)) (12) 

Near r ~ Tq eq. (6) should have this form and we obtain, for each value of ^, 
two equations: 

= E ^A^A sin(7r(//A - n))A!,{£) (13) 

A 

Tmi oc Una cos(7r(//A - T£))^a(£) (14) 
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Equations (8) and (13) can be combined to 

S\A{i))=0 (15) 

with 



S 



(16) 



[/ArAsin(7r(/XA - T^)) ; A^open ^ 
C/ata sin(7r(//A + i^Ar)) ; A^closedy 
where S stands for sine, not to be confused with the S matrix in eq. (11). 

3.2 MQDT as a Quantum Poincare map 



We now seek an alternative along the lines sketched in ref. [13] for the case of 
all channels closed. The matrix S defined in eq. (16) is the imaginary part of 
a complex unitary matrix 



^ UNAexp{m{i^A- Te)) ; A^open 



[/AfA exp(i7r(/XA + i^Af)) ; iV closed 



(17) 



Splitting real and imaginary parts in the unitarity relation 

£^£^(C'-iS')(C + iS)^i (18) 

gives 

C^C + S^S = i (19) 

C'S = S'C (20) 

Equation (20) implies that if |^(^)) is an eigenvector of S with eigenvalue 
zero, then \B{i)) = C\A{i)) is an eigenvector of 5* with eigenvalue zero, i.e. a 
solution of the transpose of eq. (15): 

S^\B{£)) = (21) 

Writing down the components of \B(i)), 

B^{£) = J2Cr,AAA{£) (22) 

A 

we see that the matrix C relates the eigenfunction for channel i in the molecular 
reference frame labeled by A to that in the laboratory reference frame labeled 
by A^. Note that eq. (22) is the same as eq. (14), except that the matrix T is 
defined for open channels only, while B is defined for all channels. 
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Next, consider the symmetric complex unitary matrix: 



= (i-25*5) +2iC*5, 



(23) 



using eqs. (19,20). This imphes that if S\A\) — 0, then 



S'S\Aj,)=I\Aj,), 



(24) 



i.e. the solutions of eq. (15) are eigenvectors of £ with eigenvalue 1. 



Similarly the I-Bat) are eigenvectors of £8'^ with eigenvalue 1 



(25) 



Eqs. (24,25) imply the determinant equation 



det(l-£*£) =0. 



(26) 



Comparing this result with the equation for the T{En) matrix defined by 
Bogomolny [9,10,11] was essential for the interpretation of MQDT as a QPM 
in ref. [13]. We will examine this relation more closely in the next sections. 

For scattering systems we have to go a little further, because we have to search 
for the phase shifts at any energy rather than for eigenvalues and eigenfunc- 
tions. For this situation we shall derive another useful form of these equations. 
Remember that in this case vectors and matrices in eqs. (24,25) depend on the 
channel number i and that the B]\f{£) are identical to the eigenvectors Tj^i of 
the S matrix. We are looking for a set of equations written for open channels 
only, whereas eq. (24) is written for all channels, open and closed. For this 
purpose we write eq. (25) in bloc diagonal form, splitting the N channels into 
open channels indexed by o and closed channels indexed by c. This is possible 
for elements indexed by the asymptotic quantum number N, e.g. S S^, not for 
those indexed by the Born Oppenheimer quantum number A, e.g. S^S. We 



obtain 



{£ £^)ooBo + {£ £*)ocBc = Bo 
{S 8^)coBo + {S S*)ccBc = Be 



(27) 
(28) 



Eliminating Be gives 



((£ S%o + {S t\c (ice - {t ^*)cc) ' Bo = Bo (29) 
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Notice then that we can rewrite eq. (17) as 

E = exp(i7rz>) JJ exp(i7r/i) (30) 

where v is the diagonal matrix in the laboratory reference frame with diagonal 
elements 

\-n ; iVopen 
vnn = < (31) 
I +t/7v ; iV closed 

and fi is the diagonal matrix in the molecular reference frame with diagonal 
elements ha- 

Then, defining 

U ^Uex.p{2infi)U^ (32) 
Equation (29) can be rewritten as 

(Uoo + Uoc (exp(-2i7r/>)ec - Ucc) ~' 4o) B„{e) = exp(2i7rr^)5„(£) (33) 

The matrix within the left large parenthesis no longer depends on Tf, it is 
defined in open channels, and for each £ has the same eigenvalue (eq. (11)) 
and the same eigenvector (eqs. (14) and (22)) as the S matrix. We thus obtain 
the expression 

S ^Uoo+ Uoc (exp(-2i7ri>)cc - Wcc) ^ Uco (34) 

for the S matrix. 

3.3 Computational advantages 

The usual MQDT (sec. 3.1) uses the equation det<S = (eq. 16) which can be 
rewritten as a polynomial in tan(7rrf) whose degree equals the number of open 
channels. Each root gives a well defined r because, according to eq. (11), the 
eigenphases of the S matrix can take values between and 1 (or -0.5 and 0.5) 
only. Depending on the choice of the total energy there are the possibilities 
that all channels are closed (there is then only one value of the set {^a}) or 
open. In the last case one obtains as obvious solutions the values = //a (see 
(eq. 16). This shows that the are exactly the usual collision phase shifts 
when the energy is high enough for all channels to be open. In the general case 
equation ( 16) is not solved by a diagonalization algorithm. The usual way to 
find a solution is 

• All channels closed: No . All matrix elements depend on energy E through 
the un- One varies E and locates the zeros of the determinant by a root 
searching algorithm. 
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• Some channels closed, some open. One fixes a given value of E and searches 
the values of for which the determinant vanishes. The search may be done 
by computing the determinant and varying like in the preceding case, or 
by writing the polynomial in tanr^ and computing its roots. In any case it 
is a root searching algorithm. 



As explained in ref. [13], where illuminating figures are given, the practical 
problem with such a standard approach is that when zeros are very close 
the determinant crosses twice or more often the zero line at extremely close 
values, and nearly tangent to the zero line. This makes it very difficult to 
avoid the missing of some pairs of zeroes. Moreover, in such a situation it is 
very difficult to recover the wavefunctions individually as they switch between 
nearly orthogonal states for a very small change of E or T£. This situation 
is quite common in the case of a nearly integrable phase space, where there 
are systematic near degeneracies associated to the lack of level repulsion due 
to approximately conserved quantum numbers. The problem of missing levels 
becomes particularly severe for the high values of angular momenta we use to 
study the semi classical limit. 



In our QPM formulation (sec. 3.2) for bound systems the relevant practical 
aspect [13] was that we could diagonalize the unitary matrix £ by standard 
techniques. Then the equation is solved by a search of zeros of the eigenphases, 
which move monotonously as a function of energy. Though it still requires a 
root searching algorithm, it is much more efficient for near degenerate lev- 
els than the standard MQDT equation det»S = and may be accelerated 
by linear interpolation. In the case of approximate double degeneracies, the 
diagonalization algorithm always reliably provides a good pair of orthogonal 
wavefunctions. 



Similarly in the ionized case the computation of the S matrix with eqs. (33,34) 
is reduced to a standard diagonalization of a complex symmetric matrix (for 
any real energy) , with an auxiliary inversion of a matrix in the closed channel 
space. The components of the resulting wave function on the open channels 
{Bo{€) = Toi) are obtained by the diagonalization and the components on 
the closed channels are then deduced with eq. (28). An other way to use this 
equation is to look for the zeros of the matrix to invert in the complex energy 
plane, which gives poles of the S matrix. We will use both techniques in the 
next sections. 
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4 Compcirison between classical and quantum evolution on the sur- 
face of section 



4-1 Principles 

One of the basic tools we use to study properties of our molecular system is to 
compare the PM with Husimi [33] or Wigner [34] plots of eigenfunctions, as 
well as the time evolution of such quantum distributions in phase space with 
that of a swarm of classical trajectories. In such plots we represent the angular 
part of the wave function, ignoring its radial part. In other words we represent 
both the classical and the quantum evolution on the surface of section, thus 
relating the PM and the QPM. 

First wc recall that we can perform the comparison of such objects in two 
different frames, the molecular or Born-Oppenheimer frame and the labora- 
tory frame. The molecular frame wave functions have coefficients labeled by 
A, the laboratory frame wave functions have coefficients labeled by N. The 
coefficients of a given wave function in those two reference frames are related 
by eq. (22). According to eq. (17) the matrix C which performs the transfor- 
mation is composed of the matrix U, which transforms between the molecular 
and the laboratory frames and of phase shifts in the two frames. We will ex- 
plain the meaning of these phase shifts at the end of this section. Note that, 
according to [35, eq. (A2.1)], the asymptotic value of U for large angular mo- 
menta reduces to a rotation of 7r/2 around OY. U thus basically performs a 
rotation between two orthogonal polar axis M and N. Wc shall focus on four 
key points in an electron orbit: perigee, apogee, and the two crossings of the 
surface of the sphere of radius ro, the border between the free motion and the 
collision region. For hyperbolic motion the apogee does not exist. Wave func- 
tions for all these four points are naturally contained in the MQDT. Perigee 
coefficients are the A\ introduced in eq. (6). They are naturally expressed in 
the laboratory frame. In the case where all channels are bound, apogee co- 
efficients are the components Bj^ of the eigenvector \B) of 5* introduced in 
eqs. (21,22). They are naturally expressed in the laboratory frame, but may 
be expressed in the molecular frame by multiplication with the matrix f/*. The 
coefficients at the crossing of the sphere of radius tq before and after collision 
are the angular parts of incoming ^p'{r) and outgoing v^^(r) in eq. (10). They 
are defined for unbound motion, but the basis for the validity of MQDT is that 
at Tq the inner part of the wave function is approximately the same for bound 
and unbound motion. So even for bound motion one can speak of "infinitely 
far from the core" motion at tq. This was a key point for the derivation of the 
classical limit of MQDT in ref. [8]. The incoming and outgoing angular parts 
are naturally defined in the laboratory frame, but inserting eq. (9) into eq. (6) 
we see that they are respectively e™i^A^ and e^^'^'^^lA in the molecular frame. 
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These considerations are clearly displayed by combining eqs. (24) and (30): 



= det(i-,f* S) 
= det(i-e''''^ C/*e' 




(35) 



As written this is the determinant of a matrix in the molecular frame. It rep- 
resents a QPM between perigee and perigee, and the associated eigenfunctions 
arc the These are the angular part of the perigee wave function. They 

depend on the energy E and (if there are open channels) on the phase shift 
Ti which cause the determinant to vanish. The second line of this equation 
reads, from right to left: apply half a collision (e"^^), rotate the axis from OZ 
to OX by il, apply half a free rotation (e*'^'^), you are at the apogee. Then 
apply half a free rotation, rotate from OX to OZ, apply half a collision and 
you are back to the perigee. 

One may cyclically interchange the matrices in the product entering the de- 
terminant. Then eigenfunctions of the associated matrix are transformed ac- 
cordingly. Moving the matrices step by step one place from left to right this 
gives first the QPM for the outgoing wave, \e™'^'^Af^ in the molecular frame, 
then the same QPM in the laboratory frame, next the QPM between apogee 
and apogee in the laboratory frame etc. This provides the explanation for the 
phase shifts inchidcd in eq. (22): going from to Bjq involves not only a 
change of the basis, but also half a collision, and half a free rotation. 

4-2 All channels closed 

That this is the correct interpretation is illustrated in Fig. 3, where we com- 
pare the time evolution of a (quantum mechanical) Husimi distribution and 
a corresponding classical swarm of trajectories for a bound system. The pa- 
rameters correspond to a mean energy for which T^/Tj^i = 1/3, so that a 
free rotation is 1/3 of a turn, and a strong coupling K = —3. 98765 tt giving a 
strong stretch around OZ, nearly ± two full turns at the ±OZ poles (and zero 
at the equator). The remaining parameters are chosen as L — AO, J — 200, 
2B = 2.5 10~^^ a.u.. The angular momenta are large as compared to those 
accessible in current molecular experiments, especially L which usually is lim- 
ited to the range 0-2 in two step laser excitations. They correspond to semi 
classical scaling of the experimental situation [6,7], as explained in section 2.3. 

The quantum time evolution is computed by first obtaining a set of 2L + 1 
consecutive cigcncncrgics Ei and the corresponding angular wave functions A\ 
for increasing principal quantum number i>j by one (or a multiple of such sets 
to increase statistical accuracy of the results) . An initial radial wave function 
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Fig. 3. Comparison between Classical and Quantal time evolution. All chan- 
nels closed. Parameters vj = 875 — 880 a. u., fi^o = 39.8765, corresponding to 
Te/TN = 1/3, K = -3.98765 TT. (a): Start, (b)-(d): Application of £. Half a free 
rotation (e^'^'^), 1/6 of a turn clockwise around OX leads from (b) to (c) (notice the 
flip from top view to bottom view), U changes axis from OX to OZ, half a collision 
(e"^^), stretch around OZ leads from (c) to (d). (d)-(f): Application of Half a 
stretch around OZ, U^, half a rotation around OX. 

is then expanded in this eigenbasis. The time evolution results by applying 
exp{iEit) to the A\. Notice that it is applied only to the angular part, not to 
the radial part of the wave functions. But for integer times, which correspond 
to a full orbit because we express energies in terms of a mean v, the radial 
part is the same, so the corresponding radial overlap is one. This procedure 
is consistent with the QPM point of view. What we neglect in doing so, is 
the spreading with time of the radial wave packet. This is consistent with 
a semi classical picture. For a full quantum computation we would need to 
choose a radial distribution /(r) peaked around the corresponding radius r, 
perigee, apogee, ro, and to develop this initial state, angular and radial part, 
in an eigenbasis. The initial quantum distribution is the Husimi distribution 
of a \LL) state. The corresponding classical distribution is a circle around OZ 
with radius given by cosOl = L/ \J L[L + 1). In Fig. 3 the top pair (a) equally 
corresponds to incoming, perigee and outgoing waves because the collision is 
a rotation around OZ and the distribution is rotationally invariant around 
this axis. The next pair (b) is the apogee. Quantum mechanically we cannot 
compute it from the previous distribution because we can compute only the 
evolution with an integral number of turns. It is thus obtained by a rotation 
of the initial distribution around OX by 1/6 of a turn. But then everything 
is computed from the four initial distributions by applying the above method. 
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One sees perfectly the evolution: 1/6 of a turn around OX (c), half stretch 
around OZ (d), half stretch along OZ (e), 1/6 of a turn around OX (f) etc. 
The quantum distribution follows perfectly the classical one in such short 
times. For longer times they begin to differ due to the finite value of h. 



4-3 First Channel Open 




Fig. 4. First channel opening. Classical PM. Since Poincare maps have a symmetry 
C2 around OX, only half the sphere is needed. We use a stereographic projection 
from +OZ of the lower part of the sphere. The open part is a small cap near the 
+OX axis. Large black squares at its border are ionizing trajectories. 

To begin with, we select an energy just above the first threshold. As for the 
rest of this section we select L = 50, J = 100, 2B = lO^^^a.u., so varies 
between 50 and 150. As compared to the bound case, the ratio L/J has been 
increased from 1/5 to 1/2 to increase the transfer of energy between electron 
and core. The PM is given in Fig. 4. The ionization region corresponds to a 
very small white cap around the +0X axis. The large black squares at its 
border are the points of ionization of some trajectories. 

The value of the phase shift r of the S matrix as a function of energy is shown 
in Fig. 5: for each energy there is only one open channel, thus one value of r. 

The Husimi plot of the wave function at the first channel opening (see Fig. 6(a)) 
is located in the open part of the sphere. Resonances for nearly bound states 
appear as increments of r by one. The widths of the resonances vary widely. 
The physical reason for this is easily displayed by Husimi plots for the reso- 
nances observed in Fig. 5. The very broad resonance which spans all of Fig. 5 
is seen to be localized in the chaotic region, overlapping widely with the open 
part of the sphere: Fig. 6(b). We see three other resonances of decreasing 
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Fig. 5. First channel opening. Phase shift. It appears at the first ionization thresh- 
old, which is negative (—872610"^°), simulating a stray space charge: see text. 
There is only one phase shift r, because there is only one channel open, in this en- 
ergy range. It increases by one for resonances between the open channel and some 
closed channels. Coupling ji^Q = 4. 

width, which are also localized in the chaotic sea, but farther from the open 
part, as shown in Fig. 6(c)-(e). Note that these resonances live on the border 
of an integrable island, and are thus of similar type as states computed in 
ref. [36]. The state 6(f) is even narrower as it hves on or around an integrable 
island within this chaotic zone. It is seen as a one point glitch in Fig. 5. In 
addition there are much narrower resonances which cannot be found easily 
by the previous method. They are localized on or beyond the large regular 
region in the middle of the sphere, and have exceedingly small width, due to 
an exponentially small overlap of the Husimi function with the open part of 
the sphere. 

Using the traditional method which looks for a zero value of a determinant we 
have been able to locate them by the following trick. We can set artificially 
a lower than zero threshold for ionization both in classical and in quantum 
mechanics: this would correspond to the existence of a long range screening 
potential, which might be due to stray electric fields or a space charge. In 
the classical simulation, this corresponds to considering an electron which 
after collision has an electronic energy higher than a given negative threshold 
^'thresh as iouized and not returning. In quantum mechanics, since the only 
difference between open and closed channels is the replacement of by — r, 
this corresponds to performing this replacement for a finite value of z/jv, which 
corresponds to this i^thresh- We have selected such a case: see the negative 
energy for opening of the first channel in Fig. 5. If we do not apply this trick, 
the whole sphere is closed for the same value of the parameters. We can thus 
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Fig. 6. First Channel Opening. Husimi plots of resonances, (a) Wave function of 
the open channel out of resonances. It is quantized on the open part of the sphere, 
(b): very large resonance which spans all Fig. 5, it has a large overlap with the 
open part of the sphere, (c)-(e) three narrow resonances seen in this figure (in order 
of decreasing widthes) at -8.7199110-^ a.u. (width 1.11 IQ-^^ a.u.), -8.7215610"^ 
a.u. (width 2.26 10"^^ a.u.) and -8.72065 10"^ a.u. (width 3.15 lO"^^ a.u.): they are 
located in the chaotic region connected to the open part, (f) the narrowest resonance 
seen in this figure, as only a one point glitch at —8.71780 10~^ a.u. (width 2.67 10~^^ 
a.u.): located on a regular island with some overlap with the chaotic sea. (g)-(h) 
extremely narrow resonances located in the regular regions: their widths are 2 10~^^ 
and 1.5 10~^^ a.u.. (i) an extremely narrow resonance located in the second chaotic 
sea between them. Its width is 1 10~^^, intermediate between the two previous ones, 
because it is separated from the open part by a regular region. 

use MQDT for closed channels to compute energies and eigenfunctions of all 
channels. We select by inspection the wave functions quantized in the regular 
regions and run the open channel MQDT with a very small step to locate the 
resonances. Thus, we have been able to locate resonances of widths 2 10~^^ and 
1.5 10~^^ a.u. located respectively in the regular regions near the equator and 
the bottom of the sphere, and also a resonance of intermediate width 1 10^^^ 
located in the second chaotic sea between these two regular regions. In any case 
we have checked that the Husimi plots for the closed level and the resonance 
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are essentially equal, which provides a good check of our interpretation. 

But this situation is quite unsatisfactory in the general case of truly posi- 
tive energy, where we cannot use this trick to locate such narrow resonances. 
Fortunately the method of locating the poles of the S matrix in the complex 
plane indicated at the end of section 3 is much more powerful. We can use 
the contour integration to locate the resonances by dichotomy in the complex 
plane, and finally use a root searching algorithm to compute them accurately 
once they have been approximately located. To check this method we have 
determined for example that the resonance at —8.71991 10~^ a.u. (Figs. 5 and 
6(c)) corresponds to a pole at complex energy value 1.11 10~^^ a.u., consistent 
with the previous method, and that the resonance at —8.71780 lO^'', seen as 
a one point glitch in Fig. 5 corresponds to a pole at complex energy value of 
2.6710-15 a.u. (Fig. 6(f)). 

4-4 Many Open Channels 

We now increase the total energy to ~ 4.50 10~^ a.u., such that slightly more 
than a half of the sphere is open. Ionization occurs in a cap around OX with 
cos a > —0.1. We select two values of the coupling, a weak /^so = 0.4 and a 
strong fjL^o — 9.81. The remaining parameters are chosen as L = 50, J = 100, 
2B = 10-1° a.u. as before. 



4-4-^ Weak coupling 

X X 




Fig. 7. Weak coupling, (a) Ordinary Poincare Map. Iteration stops when the electron 
ionizes (black squares) (b) Combined Poincare Jung map. When the electron ionizes 
iteration is continued with the Jung recipe (Fig. 1(b)). It follows circles on the open 
part of the sphere until it enters back into the closed part. 

Fig. 7(a) displays the ordinary PM. Iteration stops when the electron ionizes 
(black squares) and a new trajectory is launched in the closed part of the 



27 



sphere. Fig. 7(b) shows the combined PM+JSM. Except for the regular islands 
near the —OX axis, it was obtained by running a single trajectory, i.e. using 
the PM for negative electron energies, and the JSM for positive energies. When 
the system is in the bound chaotic part it diffuses slowly towards the back 
or the front of the sphere. This diffusion mechanism was studied carefully in 
ref. [21]. When it ionizes, the JSM generates a set of one bounce events, tracing 
a semi circle around OZ in the open region, until it enters the bound region 
at the end of the semi circle and starts anew to diffuse. Due to the small value 
of ji the number of steps needed to complete the semi circle is very large. The 
trajectory thus behaves like one injected very close to a parabolic manifold. 
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Fig. 8. Phase shifts. Weak coupling ijl^q = 0.4 

Quantum mechanically Fig. 8 displays the phase shifts for the eigenchannels. 
There are 23 channels open at this energy. We readily see that it displays a 
set of resonances. 

Some Husimi plots of the eigenfunctions of the open channels are displayed in 
Fig. 9(a-c). They follow closely the trajectories of the JSM. In a previous letter 
we noted [15], that parabolic manifolds in a JSM may lead to eigenstates of 
the S matrix localized near this manifold (see also below). As we mentioned 
above the semicircles behave quite similarly, and as we do find states local- 
ized near these semi circles this again proves the applicability of the JSM to 
the description of the ionizing part of the spectrum in a semiclassical limit. 
Fig. 9(d) displays the Husimi plot of the first resonance seen in Fig. 8. It 
is quantized along the border of the regions of bound and unbound motion, 
yielding a large overlap with the open part resulting in a large width of the 
resonance. All resonances seen in Fig. 8 have similar Husimi plots. Functions 
lying further from the open channels would again correspond to extremely 
narrow resonances, which cannot easily be seen in the eigenphases with the 



28 



Fig. 9. Husimi plot of the MQDT eigenfunctions. Weak coupling, (a-c) Open chan- 
nels are quantized along semi-circles which pertain to the Jung Map. (d) All reso- 
nances seen in Fig. 8 are quantized near the border of the separation between open 
and closed parts of the sphere. 

present resolution. 
4.4-2 Strong coupling 

Choosing as coupling strength fx^Q = 9.81 the open part slightly greater than 
half a sphere. Fig. 10(a) shows the combined PM + JSM. Except for the fixed 
points described below all points of this combined PM + JSM pertain to a 
single trajectory. Due to the large value of the coupling this trajectory jumps 
randomly in and out of the bounded part within a few steps. The combined 
map is chaotic everywhere in this part. Fig. 11 displays the phase shifts. One 
sees a lot of avoided crossings between the phase curves, signaling a chaotic 
situation. A typical Husimi plot of an arbitrary wave function is shown in 
Fig. 10(b). 

The important fixed points of this map (apart from isolated unstable fixed 
points embedded in the chaotic sea) are 

(i) the —OX axis in the bound part (PM) surrounded by a very small island 
of stability, too small with the present value of h to capture a single 

quantum state. 

(ii) the ±0Z axis in the open part surrounded by a larger regular island of 
the JSM. These are fixed points because the collision induces a rotation 
around OZ, entirely embedded in the open part (JSM). 
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Fig. 10. Combination of PM and JSM. Strong coupling. Left: top view, right: 
bottom view. We do not use the stereographic projection as in Figs. 4, 7 because 
the parabolic manifold is badly seen with it. (a) Classical map. Small dots: PM, 
larger dots: JSM. The open part is slightly larger than half a sphere. The PM has 
only a very small island of stability at the back, too small to capture a quantum 
state with this value of h. The JSM has two regular islands around ±OZ, and the 
parabolic manifold of fixed points on the equatorial plane. All the remaining points 
lie in a single chaotic sea obtained with only one starting point, iterated by the PM 
when bound, and by the JSM when ionized, (b) Husimi plot of a typical arbitrary 
state: it extends randomly over the whole sphere (c) Husimi plot of a quantum state 
quantized on the parabolic manifold, (d) Husimi plot of a quantum state quantized 
on the zizOZ cap. 
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Fig. 11. Eigenphases as a function of energy. Strong coupling /iso = 9.81 

(iii) all points in the equatorial plane in the open part (JSM). An electron 
entering the collision region with 6^ = vr/2 and arbitrary ifi, up to the 
restriction that the electron energy should be positive, will emerge from 
the collision region with the same angular momentum direction because 
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the increment SipL = Kcos9l equals zero. We thus have a parabolic 
manifold of fixed points, 
(iv) if \K\ is chosen larger than n ■ 2tt, that is /x^, > ^ we would also find 
lines mapped onto themselves at angles 9l corresponding to KcosOl — 
0, ±27r, • • • , ±n • 27r 

In particular the parabohc manifold is embedded in the chaotic sea. A nearly 
horizontal curve may be observed in the middle of Fig. 11. It corresponds to 
the fixed points of this parabolic manifold in the equatorial plane as shown 
by Fig. 10(c) (see also ref. [15]). For e.g. the stadium billiard, such manifolds 
of fixed points embedded in a chaotic sea have significant quantum effects. 
There they give rise to the so called bouncing ball states which in turn lead 
to deviations of the spectral fluctuation properties in the eigenvalue spectrum 
from the random matrix results [37]. We expect similar quantum effects of 
the parabohc manifolds on statistical properties of the eigenphases of the 
scattering matrix. When considering eigenphases of the scattering matrix as 
a function of the energy E or the coupling strength yU^, we flnd eigenphases 
which merely change. The other fixed point near the ±0Z poles, surrounded 
by a significant island of stability in the JSM, is less conspicuous in Fig. 11. 
It corresponds to a horizontal line near r = —0.3, barely visible due to strong 
avoided crossings. It nevertheless traps a quantum state as seen in Fig. 10(d). 



5 Comparison between semi-clcissical and quantum evolution on 
the surface of section 

5. 1 Principle 

We want to compare our results with the original work of Bogomolny [9,10,11], 
which was based on the semi-classical formula 

T{^, ^ ) = E V exp ^-S{^, ^ ) j , (36) 

where, as in the remainder of this section, we use the notations of the review 
paper of W.H. Miller [38]. In particular the i and the signs in the argument of 
the square root function take care of the Maslov index. 

We first notice that the simple minded way of using this equation, computing 
S = J pdq using the formulae obtained in the previous classical part of this 
paper and the previous ones does not work, because our "classical" molecular 
reference frame is a moving frame, implying that Lz^ is not conjugate to <^£^. 
Indeed consider first the collision step. The change of orientation of L in this 
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frame entails, by conservation of angTilar momentum, a change of the orien- 
tation of N, and thus of the orientation of this frame, by an angle called 5ip']^ 
in section 2 and ip% in Appendix A.l. This extra rotation of the frame im- 
plies that Lzc is no longer the generator of rotations of around the OZc 
axis. For the collision step this problem can be solved by using the "quan- 
tum" molecular reference frame, which is fixed, since the collision is taken 
as instantaneous, and that this frame is obtained from the laboratory frame 
by a rotation with fixed Euler angles <^m,6'm,0 (see Appendix A). That this 
analysis is correct is shown by noting that the transformation produced by the 
collision in the "quantum" reference frame t^f/j^ = KcosOl is canonical (the 
exterior differential product dLzgAdip^ is conserved), while the corresponding 
for the "classical" frame quantities is not. Moreover the transformation from 
the "quantum" to the "classical" molecular reference frame is not canonical as 
shown in Appendix A.l, where we have computed the momentum conjugate to 
(fil, such that this transformation becomes canonical (eq. A. 7). Nevertheless, 
since our goal is physical evidence we prefer to stick to our previous practice of 
using Lzc? V^L- The only drawback of using non canonical coordinates is that 
the map is not area preserving on the sphere, a minor effect for L/J=l/5or 
1/2 as chosen in this paper. 

For the collision step we thus can solve this problem by using the "quantum" 
molecular frame, but the problem reappears for the Coulomb step, where this 
frame is also moving. We could compute everything in the laboratory frame, 
but it is more meaningful to use the "quantum" molecular frame for the col- 
lision step, the laboratory frame for the Coulomb step, and two semiclassical 
formulae for the change between these coordinate frames. This procedure fol- 
lows closely the purely quantum formulation, giving a semiclassical approx- 
imation for each of the steps which enter in the discussion of section 3, as 
shown in Table 1. It provides direct evidence that MQDT indeed is the exact 
quantization of the Bogomolny formulation, in particular that eqs. (26, 35) are 
the exact counterpart of the semi classical eq. (1.20) of ref. [10], remembering 
that £ describes half a turn, while T describes a full turn. 

As Poincarc section we select in coordinate space the intersection at radius 
tq with outgoing trajectories, i.e. an "after collision" Poincare map. We thus 
write: 

T ^VU^^U^U^^ (37) 

where from right to left the C/"s indicate the unitary semiclassical operators 
for the transformation from the "quantum" molecular reference frame to the 
laboratory frame, free motion in the laboratory frame, back to the "quantum" 
molecular frame and collision in the "quantum" frame. All these operators are 
built with F2 I F3 kind of generating functions, i.e. they are of the form 
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d^F:,{I'.^-)/Drd'^ 



exp ^^{Fs{I',ip)) 




(38) 



where /' is the momentum conjugate to the coordinate </?' and the resulting 

integral is evaluated with the stationary phase approximation. This generat- 
ing function is of type as it is a function of the old momentum /' and 
the new coordinate (p. An auxiliary Fourier transform is used to go from old 
coordinate if/ to old momentum /'. Evaluating the inner d/' integral by the 
stationary phase approximation amounts to computing an Fi type generating 
function, and inverting the order of the Fourier transform and of the F integral 
amounts to defining an F2 generating function. Note that we have defined by 
this double integral a U operator which goes from coordinate (/?' to coordinate 
(fi representation, while the F3 integral on its own gives an operator from the 
old momentum I' to the new coordinate (f. Hence, we can compute eq. (37) by 
performing four integrations instead of eight by alternating between the co- 
ordinate and the momentum representations. We can compute with the same 
procedure a momentum coordinate Taa' matrix by simply introducing two 
Fourier Transforms to the left and the right of the operator T. As momenta 
are quantized Tw depends on discrete variables, while T{ip, ip') is a function 
of continuous arguments If there are no singularities all these possibilities give 
the same result. We shall use them freely in order to obtain generating func- 
tions which we can compute analytically and we shall see in section 5.3 that 
this choice offers possibilities to circumvent singularity problems. 

5.2 Computations 

First notice that the conjugate pair of coordinates in the "quantum" molecular 
frame is the known {L^^,!^^} pair. The momentum in the laboratory frame 
is the modulus N of N. Intuitively the conjugate angle is the angle <^jv of 
rotation of the core axis M around N, referred to a reference axis which only 
can be J. Explicit formulae for computing the transformation are given in 
Appendix A. 2, where we also show that the transformation from {Lzq, V?^} to 
{iV, 93jv} is indeed canonical, justifying our intuition concerning the definition 
of lpn- 

The computation of the F3 generating functions for the collision and the free 

rotation steps follows the same scheme. In both cases the momentum is con- 
served, while the angle is increased by an amount which is a function of the 
corresponding momentum / only 5^9 = /(/), respectively ^(^^ = ^^Zq and 
5ipN — —SP — -|-27rTe/Tjv (the free rotation is direct in the laboratory frame 
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and seen with inverse orientation in the "classical" molecular frame), which is 
a function of only for a given total energy E. Using the equations defining 
F^ll'jip), namely ip' = dF^/dl' and / = dF^/dip we obtain 

Fs{I',ip)^cpI'- j f{l')dl' (39) 

i.e. using [8, eqs. (A11,A12)]: 



2L 



(40) 



The computation of the F2 or F^{Lzqi^n) is more tricky. First, whether it is 
a F2 or F3 generating function depends on which is the original and the new 
frame, and we need both transformations. To go in the opposite direction we 
use the unitarity of the operator U . We start with the following formula which 
expresses N as a function of the independent variables: 



^J^- {Lzq/ sin cpNy±^L^ - {Lzq/ smifN)' 



(41) 



This formula can be computed analytically by eliminating between eqs. (A. 11) 
and (A. 15) and using eq. (A. 13). It is more satisfying to notice that a little 
geometrical thinking on Fig. 2(b) shows that it results from projecting the J 
and L sides of the triangle given by the relation J = L + N onto the N side, 
because the common height is equal to \Lzq/ sin 997V |- By integration over (^jv 
we obtain 

HLzq, ^n) = {Lz^AtcXY (^J2-(L^^/sin^^)2, Lz^ cot if^) (42) 
+J ArcXY (^J cos ipM, \J J'^ ~ {-^Zq/ sin(^jv)^sin(^jv)) 
± (Lzq ArcXY (yZ^^^^L^^Tsin^ivj^, Lzq cot ipn) 
+L ArcXY ^Lcosf^AT, y^L^ — {Lzq/ siny^Tv)^ siny^Ar)^ 

where ArcXY is an unambiguous notation for the angle of a vector whose 
projections on the X- and F-axes are known (the order of parameters for the 
analogous atan2 used in some programming languages is not fixed). 

Notice that the derivative of F2 with respect to Lzq, which must give y?^, 
equals the sum of the coefficients of Lzq in the first and third terms of eq. (42) 
(this also has a geometrical interpretation). This shows that the sum of these 
terms is precisely equal to 'pj^Lza- Using this result we can solve the problem 
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of the integration constant depending on Lzq , which appears in the integration 
over (/9jv. We follow continuously the various angles over ranges greater than 27r 
during the computations in order to compute the exact values of the various 
ArcXY. 

Next wc compute the actions S{'^, '^') and S'(A, A'). We start from a pair of 
conjugate {L'^^,ip^ }, i.e. a point on the Poincare sphere, and compute ana- 
lytically the final Lzq, (fl as a function of these. To compute, say S{(p1, ipf ), 

we first invert numerically the relation between ip^ and L'^ for fixed ipl . 
Usually, there are several solutions. For each solution we then compute S* as a 
function of (f^ and L'^^ by adding the various F„ generating functions which 
enter in the exponentials of the U operators which are used to build T. 

Wc do not need to evaluate all intermediate integrals by stationary phase to 
obtain the amplitudes. Indeed, it was shown in ref. [38] that the formula for 
the amplitude given in eq.(36) is a direct consequence of the action being a Fi 
generating function of a canonical transformation between the initial and the 
final states. We thus only need to compute the second derivative of S with re- 
spect to (f'^ and ifl ■ For this we use the derivatives of S with respect to initial 
variables L'^^ and , that we have computed analytically by following the 
various steps, and general calculus formulae for second derivatives of inverse 
functions. We have checked that these amplitudes obey the relationship. 



(43) 



which was shown in ref. [38, eq. (2.54b)] to be a simple consequence of the 
conservation of probability. This is a very useful test of our way to compute 
S since the second term of this equation does not depend on S. 

5.3 Singularities in the semi classical computations 

We then use the resulting T(ip, cp/) and ^(A, A') to compute the evolution of 
wave packets as was done in section 4.2. The first problem which appears is 
the existence of divergences of the amplitudes. Fig. 12 plots the amplitudes 
of T as a function of the initial position on the Poincare sphere, which was 
used as a first step to compute them. The black fines of maxima are actually 
divergences. They are at different positions for the two representations. 

The physical meaning of these divergences is the existence of focus singulari- 
ties. We discuss the (p representation first, by referring to Fig. 13(a). Various 
trajectories starting from the same coordinate point with different momenta 
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Fig. 12. Divergence of the semiclassical amplitudes, (a) T{ip, ip'). (b) T(A, A'). The 
part of a straight Hne going through the OZ pole in the lower left of the sphere is 
the input corresponding to Fig. 13(b): the fact that it is slightly distorted is due 
to plotting in the "classical" reference frame instead of the "quantum" reference 
frame. 

refocus at the same point (on the map). On the sphere the coordinate is </?' 
and the momentum is L^. The points aligned on the straight line in the lower 
left part of Fig. 12(a) and (b), have the same value of ip' and a varying L'^. 
They thus constitute a fan of initial trajectories. The image of this line is given 
in Fig. 13(b). The focus is caused by the folding part of this image, precisely 
the point where it is tangent to a radius from the OZ pole: same position 
(f independent to first order of the initial momentum L'^. The divergence of 
the amplitude is a consequence of the presence of the term d(f/dL'z\(fi' in the 
denominator of the second part of eq. (43). 



X 




Fig. 13. Focus singularities, (a) represents such a singularity in an ordinary (flat) 
coordinate space, (b) plots the image of the straight oblique line in the lower left 
part of Fig. 12 (a) and (b) which cross the singularity line in the T((p, p') representa- 
tion. All starting points on this segment have the same coordinate ip' , but different 
momenta L^, i.e. they are equivalent to the starting points in (a). The back folding 
of the image in (b), common especially to chaotic situations, is the cause of the 
focus singularity: see text. 

The problem of singularities is common in any semi classical analysis. Their 
solution, usually called "regularization" follow two main routes. The first, com- 
mon in physics textbooks, is to revert to a full quantum computation near the 
singular point, and to match the results with the semi classical computation 
distant from it. Notice that the MQDT analysis is such a fully quantum so- 
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lution. The second, common in some mathematical works, is to notice that 
we can circumvent the problem by changing from coordinate to momentum 
representation through a Fourier transform. The singularity is a result of a 
projection from phase space onto coordinate space; there are no singularties 
in phase space itself. The back bending in Fig. 13(b) which is tangent to a 
meridian of the sphere is not tangent, at this point, to a parallel. The starting 
points in Fig. 12 which project onto the same ip' project to different L'^. The 
lines of singularities exist in both representations in Fig. 12, but they appear 
at different points. So one can compute the points of the map by different 
routes, a possibility rooted on the possible choice at each U step in the com- 
putation of T (eq. (37)) to select the tp or Lz representation, both in the initial 
and the final state. The only remaining problem is the fact that both lines of 
singularities share the ±0Z poles. To solve it we would have to use a more 
general change of canonical coordinates on the sphere, a rotation which can 
change the position of the poles. We do not attempt to do this in the present 
paper. 



Start Class. MQDT T{ip,(p') T{A,A') 




Fig. 14. Effect of focus singularities. The starting points lie on a circle. We show 
the first iteration by the Poincare map, and the corresponding Husimi plots with 
different techniques. First row: starting in a point without singularity. Second row: 
starting in a point with a (p focus singularity, at the intersection of the circle with 
the straight line in the lower left part of Fig. 12(a) and (b) and the singularity line 
in the (p representation. 

Illustration of this is given in Fig. 14. The starting point in the first row is 
selected to correspond to no singularity in any of the two representations. All 
computations, classical, quantum and semi classical in both representations 
agree nicely. In the second row the starting point lies on the singular line in the 
p representation, but on a regular point in the A representation. The result 
of the T{p,pf) computation is highly distorted, while that of the T(A, A') 
computation agree nicely with classical and MQDT computations. 
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6 Conclusion 



We have formulated the multichannel quantum defect theory for Rydberg 
molecules as the quantization of a classical symplectic map composed of a 

standard surface of section for negative electron energies and a Jung scatter- 
ing map for positive ones. This formulation elucidates MQDT from an entirely 
different point of view. Simultaneously we have shown it to provide consider- 
able computational advantages. 

The method was applied in the context of a simplified model, where only one 
rotational degree of freedom for the molecule was taken into account. Simulta- 
neously, the absolute value of the electron angular momentum was fixed. This 
approximation allows for a two dimensional surface of section, and simplified 
calculations, but is not essential to the representation of MQDT as a quantum 
Poincare map. These applications demonstrate the simplifications achieved. 
Finally we establish the relation to Bogomolny's [9,10,11] semi classical for- 
mulation of the quantum Poincare map. Considering caustics we can clearly 
see the basic difference between the exact and the semi-classical map. 



A Reference Frames 

We use three reference frames, the laboratory reference frame Oxyz, and 
two molecular reference frames, called "Classical" OXcYcZc and "Quantum" 
OXqYqZq molecular reference frames. 

In purely classical computations the "Classical" right handed molecular refer- 
ence frame is defined in such a way that OZq lies along the inter nuclear axis 
M of the molecular ionic core, and OXq lies along its angular momentum N. 
This is adequate because the two steps of the dynamics, namely the collision 
and the free rotation, are rotations about theses itwo axes. We always use it to 
display our results, as the physical content is most obvious within this frame. 

However, this is not the molecular reference frame selected in usual definitions 
of quantum wave functions. There the molecular frame OXqYqZq is defined 
as the one deduced from the laboratory frame Oxyz by a rotation of three 
Euler angles (with the conventions of [18]). The first two are well defined, 
namely ipM, Om, the polar angles of M. The third is arbitrary for a diatomic 
molecule, since it corresponds to a rotation around the molecular axis. It is 
given an arbitrary value or 7r/2, depending on authors; we will select as in 
[2]. 

Both reference frames share the same OZ — M axis, but have different OX 
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axes. The classical OXc — N axis is not well defined in quantum mechanics, 
since N is a momentum conjugate to the angles (fM, and 6m, which are the 
variables defining the rotational molecular wave function in eqs. (2-3). Thus 
the classical reference frame defines the usual quantum wave functions only 
in the semi classical limit. But conversely, the quantum reference frame is 
also well defined in classical mechanics. So, if we want to compare classical 
and quantum results, we do the quantum calculations in the usual quantum 
reference frame, and then plot the resulting Wigner or Husimi distributions 
in the classical reference frame, by applying the well defined classical rotation 
which brings the quantum reference frame onto the classical one. 

We shall now compute the formulae which connect these three frames. They 
are based on the fact that we can compute the position of J and N in 
OXqYqZq once we know the position of L in this frame, using simple geomet- 
rical relations: see Fig. 2(b). We first need to define precisely the laboratory 
reference frame. Since everything is invariant with respect to a rotation of 
the whole space, we may choose the Oz axis along J (Fig. 2(a)). Quantum 
mechanically this amounts to studying the | J, Mj) = \J, J) state. Conversely, 
once we know the position of J in the molecular frame, we know the position 
of the Oz axis of the laboratory frame in OXqYqZq, but there is no way 
to locate the Ox and the Oy in this molecular frame. This is a consequence 
of the rotational invariance of the problem. The corresponding angle (pM is 
conjugate to J^, and is one of the variables we have eliminated when reducing 
the dimension of phase space from 10 to 4 (or 6 if L is not kept constant). 

The position of J in OXqYqZq is a simple consequence of our choice of as 
the third Euler angle of the rotation which brings the laboratory frame Oxyz 
onto the "Quantum" molecular reference frame OXqYqZq. These three Euler 
angles are lpm, 6m, and 0. The first Euler rotation of the angle keeps J = 
Oz fixed and brings Oy onto OYq. The second Euler rotation around OYq by 
an angle brings Oz onto OZq. Conversely the position of J in the molecular 
frame is obtained by a rotation of — around OYq. It is thus located in the 
OXqZq plane with coordinates Jxq = —J sin 6m, Jyq = 0,Jzq — Jcos6m, 
always directed towards negative OXq since < 6'm < tt. 



A.l Transformations between "Classical" and "Quantum" molecular frames 

Let us compute N in OXqYqZq. For a diatomic molecule, the momentum 
N is perpendicular to the inter nuclear axis M, since the moment of inertia 
around the inter nuclear axis is zero. Thus N is in the OXqYq plane. We then 
write N = J — L in this frame 
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Nxq = — J sin 6m —L sin cos </?f 

Nyq^ -Lsin^f sin(^g (A.l) 



Nzq — JcosOm — I/Cos( 



^Zn — leads to 



cos^M^^cos^f (A.2) 



sin 9m = - cos^ 9m, 



because < 6*^^? < tt. This gives the position of J once L is known in 
OXqYqZq. Then N hes in the plane OXqYq at an angle 



^% = ArcXY {Nx^,Ny^) (A.3) 
= ArcXY (^-'\/j2-L2cos2 9^ - L sin cos ip2, -L sin sin iff 



where we have used eqs. (A.l) and (A.2), and where ArcXY is defined after 
eq. (42). 

Thus the "quantum" molecular reference frame angles for L are related to 
their corresponding "classical" reference frame value via: 



0^ = 01 (A.4) 

C Q Q 



The reverse is not obtained by simply inverting this relation. We rather express 
the relation J = N + L in the "classical" reference frame, where we know that 
N is along OXq. 



Jxc ^N+L sin cos ipY^ 
Jyc — +L sin 9]^ sin (^^ 

Jzc^ +Lcos9l (A.5) 



Since the OZ axis is common in the "quantum" and "classical" molecular 
reference frames, relations A.2 are conserved with 9^ replaced by 9^. Then 
we recall that in the OXqYqZq reference frame the projection of J into the 
OXqYq plane, which is identical with the OXcYc plane, is always directed 
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along the negative OXq axis, which locates the OXq axis in the "classical" 
frame. Thus 

"Pl = ^L- ArcXY (-N - L sin cos -L sin sin ^pf) (A.6) 
= - ArcXY - L2(l - sin^ 9^ cos^ -L sin sin 

where the second line is obtained by using eq. (1) which was written in the 
"Classical" molecular frame. 

Finally notice that the transformation between { Lzq,^^ } and { Lz^^'P^ } 
is not canonical, as can be checked by mere differentiation from the non con- 
servation of dLz Adi^L (where A denotes the exterior differential product). We 
have computed a generating function, and it can be checked by differentiation 
that a momentum conjugated to would be 

I{(p1,9l) = Lcos^L+ 

{.P L^-)F (^-0,. - sin^ ,f )E (f 

cos ^^\J J'^ - LP' sin^ 

(A.7) 

where E and F are elliptic functions, instead of = L'z = Lcos9l- It can 
be checked that the extra term vanishes when L/ J <^ 1, as it should do since 
in this case the frame recoil is negligible. The effect is thus small in this paper 
where we have selected L/J = 1/5 or 1/2. But this formula is too unwieldy 
to be of any use, especially since our main goal in using "classical" reference 
frame coordinates is physical visualizing of the results. 

A. 2 Transformation between "Quantum" molecular reference frame and lab- 
oratory reference frame 

The canonical momentum and angle in the "quantum" molecular reference 
frame are Lzq and ipl- The canonical momentum and angle in the laboratory 
frame (for the reduced phase space) are N and its conjugate angle (fN, to 
be defined, and not to be confused with ip^ (eq. A. 3), which is the angle 
of N in the plane OXqYq. We first express as a function of "quantum" 
molecular coordinates (eq. (1) was with "classical" molecular coordinates). 
Using eqs. (A.l, A. 2) we obtain 

N^^J^ + L'- 2L|^ + 2 ^L^ - Ll J J^-Ll^ cos cpf (A.8) 
Intuitively (we will justify it below) the angle conjugate to N is the angle ipN 
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of rotation of the molecular axis M around N. The reference for this angle 
can only be the position of J. Notice that the relative positions of J and N 
have nothing special. We choose as reference axis around N the unit normal 
n to the plane {N, J}: 

N X J L X J , . ^, 

" = li?M = Mj| 

(using N = J — L, X denotes the skew vector product) and the unit binormal 

m = - X n (A.IO) 

so that, after a httle algebra 



(^jv = ArcXY(n- M,m-M) (A.ll) 
^(LxJ)z^ -NJzA 



ArcXY 



|L X J| ' |L X J| 



with 



-JL^ ~ L^JJ^ ~ L2 sin -NLy 

: ArcXY I ^ y ^, I (A. 12) 

' L X J L X J ' 



That this is the correct definition of 99 at is checked by showing by mere deriva- 
tion that dLzQ A dc/?^ = 

Reciprocal relations are obtained from the right part of eq. (A.ll) 



and by extracting cos^?^ from eq. (A. 8) and sin<^^ from the left part of 
eq. (A.ll), and replacing Lzq by the previous expression, which leads to: 



(^2 = ArcXY + ^^^|L X i\\ - cos(^^|L x J|] (A.15) 
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